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SIGNIFICANCE  AND  EXPLANATION 


Spaces  of  piecewise  polynomial  functions  on  regular  partitions  have 
important  applications  to  finite  elements,  surface  fitting,  etc.  However,  the 
structure  of  such  spaces  is  rather  complicated  when  some  global  smoothness  is 
demanded  of  the  pp  functions.  For  example,  the  order  of  approximation  to 
smooth  functions  is  not  known  even  when  the  partition  is  very  simple. 

Moreover,  so  far  there  is  no  general  criterion  for  deciding  when  such 
spaces  can  be  used  to  approximate  arbitrary  functions.  We  have  conjectured  in 
[BDH]  that  a  necessary  and  sufficient  condition  for  this  is  that  the  space 
contain  a  stable  and  local  partition  of  unity.  In  this  work,  we  prove  this 
result  for  a  general  class  of  univariate  function  spaces  which  abstract  the 
properties  of  pp  function  spaces.  This  should  set  the  world  of  approximation 
theory  on  its  ears. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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,  ABSTRACT 

We  show  ^that  for  certain  translation  invariant  spaces  S  ,  a  necessary 
and  sufficient  condition  for  the  eventual  denseness  of  the  corresponding 
scaled  spaces  S^  is  that  S  contain  a  stable  and  locally  supported 
partition  of  unity. 

These  results  have  been  motivated  by  recent  work  on  approximation  by 
multivariate  piecewise  polynomials  on  regular  meshes. 
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PARTITIONS  OF  UNITY  AND  APPROXIMATION 
C.  d*  Boor1  and  R.  DeVore1'2 

1.  Introduction,  we  are  interested  in  the  existence  of  partitions  of  unity  in  certain 
tranalation  invariant  apacaa.  Ha  are  motivated  by  recant  investigations  of  approximation 
by  multivariate  piecawiae  polynomial  functiona  on  regular  partitiona.  For  example,  in  [BD] , 
two  partitiona  ware  singled  out.  The  firat  ia  the  partition  t  of  R2  into  aquarea  of 
unit  aide  length  whoee  verticea  are  the  integer  lattice  points.  The  other  is  the 
triangulation  of  S2  obtained  by  Inserting  the  northeaat  diagonal  into  each  aquare.  These 
partitions  have  acquired  the  names  "two-direction"  and  "three-direction”  mesh, 
respectively,  since  the  lines  generating  them  have  two,  respectively,  three  distinct 
directions.  The  "four-direction"  or  "criaa-croaa*  mesh,  obtained  from  the  "three-direction" 
mesh  by  drawing  in  the  other  diagonals  has  also  been  studied  [DM] • 

It  has  become  clear  from  these  inveatigationa  that  the  description  of  the 
approximation  properties  of  pp  functions,  even  on  these  simple  meshes,  is  much  more 
complicated  than  in  the  univariate  case. 

We  shall  say  that  a  partition  r  of  K*  is  regular  if  x  ia  invariant  under 
translation  by  any  a  €  Zs  ;  l.e.,  t+a«  T  for  such  a  .  We  denote  by  n£( t)  the 

space  of  pp  functions  f  with  subdivision  T  of  total  degree  k  and  global  smoothness 
p  ,  that  ia  f  €  C°(RS )  .  A  problem  of  particular  Interest  has  been  to  determine  the 
approximation  order  to  smooth  functions  (say  C*  )  from  the  spaces  s£(hT)  •  So  far  this 
problem  has  been  completely  settled  only  when  T  •  Z  as  above  [BD]  or  for  certain  special 
choices  of  k  and  p  for  the  three-  or  four-direction  mesh  [BD] ,  [BH ^ _2 ] .  [DM],  [j] .  These 
investigations  show  that  the  approximation  order  depends  in  a  subtle  way  on  the 
smoothness  p  and  the  number  of  distinct  directions.  For  example,  in  contrast  to  the 
univariate  case,  an  increase  in  p  usually  decreases  the  approximation  order  (and  not  just 
the  coefficient);  also,  an  increase  in  the  number  of  distinct  directions  may  increase  the 
order. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1 . 

2This  material  is  based  upon  work  supported  in  part  by  the  National  Science  Foundation 
under  Grant  No.  DMS-8320562. 


Lower  bounds  for  the  approximation  order  are  usually  obtained  with  the  aid  of  explicit 
methods  of  approximation,  known  as  quasi-intexpolants.  Tnese  use  special  splines  with  local 
support  known  as  box  splines.  It  is  not  necessary  to  go  into  the  properties  of  these 
splines  here.  We  only  mention  the  fact  that  it  is  easy  to  rake  up  a  partition  of  unity 
consisting  of  box  splines  and  this  is  the  important  property  used  in  the  construction  of 
guasi-interpolants. 

It  has  not  been  clear  to  us  whether  this  represents  the  only  viable  approach  to  the 
construction  of  suitable  approximants .  For  example,  in  [BOH],  we  have  asked  the  following: 

(1.1)  Is  U^o  *£(hi)  dense  in  CgfK8)  if  and  only  if  *£(t)  contains  a  good 
partition  of  unity? 

By  a  good  partition  of  unity  we  mean  a  (countable)  collection  ♦  of  functions  with 
the  properties  that 

(i)  sup^j  diam  supp  <p  <  <»  ; 

(1.2)  (ii)  <P  -  1 

(iii)  |<p|  <  -  . 

Actually,  it  is  easy  to  see  that  the  existence  of  a  good  partition  of  unity  is 
sufficient  for  denseness.  Namely,  if  x^  €  supp  <p  ,  then 

gives  an  operator  L  whose  dilate  Lj,  ,  with 

(Ohf)(x)  :=  f (hx)  , 

provides  an  approximation  I^f  in  n£(ht)  which  converges  to  f  as  h  — >  0  . 

It  turns  out,  but  is  much  more  difficult  to  prove,  that  the  existence  of  a  good 
partition  of  unity  is  necessary  for  denseness.  Here  we  shall  present  a  proof  of  this  fact 
for  the  univariate  case  but  for  much  more  general  spaces  S  than  just  pp  ones,  in  an 
attempt  to  extract  the  essential  features  of  pp  function  spaces  on  regular  partitions. 

The  ideas  of  the  proof  seem  to  carry  over  to  the  multivariate  case,  but  the 
constructions  involved  are  much  more  involved.  For  this  reason,  we  feel  that  it  is 
beneficial  to  present  here  the  basic  ideas  in  their  simplest  setting  and  report  elsewhere 
on  their  multivariate  analogues. 
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2.  The  univariate  case.  We  let  S  be  a  space  of  functions  on  K  with  the 
properties: 

(1)  S  is  translation  invariant,  i.e.,  f  €  S  "■■>  f(*±1)  €  S  ; 

(ii)  u  :»  din  S |  r n  1 1  <  “  > 

(2.1) 

(iii)  dim  S  «  -  i 

(iv)  S  is  closed  under  uniform  convergence  on  compact  sets, 
we  note  for  later  reference  that  (i)  and  (iii)  imply 

(v)  dim  S|  [0(OB)  -  »  . 

These  are  the  essential  properties  of  pp  spaces  over  regular  meshes.  However,  it  is 
easy  to  give  examples  of  such  spaces  S  which  are  not  pp.  For  example,  if  ♦  is  a  finite 
collection  of  compactly  supported  functions,  then  the  closure  under  uniform  convergence  on 
compact  sets  of  the  span  of  the  translates  :“  <p(  *  -  j )  ,  j  €  Z  ,  ip  f  ♦  ,  has  these 

properties.  The  results  of  this  section  will  show  that  any  space  S  satisfying  (2.1)  can 
be  obtained  this  way  modulo  some  finite  dimensional  space. 

Associated  with  S  we  have  the  dilated  spaces  Sj,  :“  °t/hS  ' 

Theorem  1.  If  Sh  is  dense  in  Cg(Z)  ,  then  S  contains  a  good  partition  of 

unity. 

In  preparation  for  the  proof,  we  now  construct  certain  compactly  supported  functions 
in  S  .  Let  s”  denote  the  subspace  of  S  consisting  of  those  functions  which  vanish 
on  (-•,0J.  It  will  follow  later  that  S”  is  nontrivial.  For  now,  we  prove: 

Lemma  2.  For  each  Integer  n  >  0  ,  there  is  an  such  that  any  <p  €  S  vanishing 

on  I-Rn,0)  agrees  on  [0,n]  with  some  function  from  S~  • 

Proof.  For  r  ,  j  >  0  ,  we  let  Lr>j  denote  the  space  of  functions  defined  on 

(0,j)  which  are  restrictions  of  functions  in  S  which  vanish  on  l-r,0]  .  Then  L_  .  C 

r,  3  — 

Lft  j  whenever  r  >  r1  .  Hence  (dim  Lrj)".^  is  a  decreasing  sequence  of  (nonnegative) 


integers,  therefore  is  eventually  constant.  This  implies  that  we  can  find  an  increasing 


for  all  r  >  R 


3  ' 


3-1,2, 


sequence  (Rj)  such  that 

Lr,j  “  Lj  Hlj.j 
Now  given  n  >  0  and  if  f  Ln  ,  we  define 

n^^  n  +  i  ,  rA  max  {R^n^}  ,  i»0,1,  ...  . 

Then,  since  r1  >  R^  *  R„  ,  there  is  an  f1  €  S  which  vanishes  on  [-r  ^ ,  0]  and  agrees 
with  ip  on  [0,n]  .  Similarly,  r2  >  Rn^  >  R^  and  so  there  is  an  f2  S  S  which  vanishes 
on  [-r2,0]  and  agrees  with  f.,  on  [O.n-,]  .  In  this  manner,  we  obtain  a  sequence  (fj) 
such  that  fj  vanishes  on  [— r* ^ , 0 3  and  agrees  with  f^  on  [0,n^]  for  all  i  <  j  . 

Hence  f j  converges  uniformly  on  compact  sets  to  some  f  ,  necessarily  in  S  because  of 
(2.1.iv),  with  the  desired  properties.  ||| 

Now  let 


be  the  constant  of  Lemma  2.  We  improve  Lemma  2  in  the  following  way: 


- - 3.  If  f  €  S  vanishes  on  t-R.O]  ,  then  ip  :«  f X [o , ~)  18  ln  s  ’ 

Proof.  From  Lemma  2,  we  know  an  fg  €  S*  for  which  f  -  fg  vanishes  on  [-R,1]  . 
Hence  there  is  some  fj  with  f  -j  ( «+1 )  €  S-  ,  therefore  also  f(  f  S  ,  for  which  f  -  fg  - 
f1  vanishes  on  t-R,2]  .  In  this  manner,  we  construct  a  sequence  (f  j )  in  S  for 
which  f g  *  ...  *  f j  converges  to  ip  .  Since  S*  is  closed  under  uniform  convergence  on 
compact  sets  (by  (2.1.iv)),  tp  €  S'  .  ||| 

Next,  we  construct  a  basis  of  compactly  supported  functions  for  S~  .  Recall  the 
definition  u  :=  dim  s|(0,1]  and  se^ 

M  :=  yR  . 

T  —a  4.  If_  f  €  S”  ,  then  there  is  a  function  if  E  S  which  is  supported  on  10, M] 
and  agrees  with  f  on  (0,1)  . 

Proof.  Set  J  : “  [M.M+R]  and  consider  the  functions  fj  f('-j)  ,  j=0,...,M  . 

Since  dim  Sij  <  M  ,  these  functions  are  linearly  dependent  on  J  .  Hence  '■.here  are 
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not  all  zero,  such  that 


numbers  c  ^  , 


*0  !"  1  °j  fj 

vanishes  on  J  .  Let  1  be  the  smallest  value  of  j 
ip0<*-i)/c^  agrees  with  f  on  [0,1]  and  vanishes  on 
Therefore,  by  Lemma  3,  i|^  :»  *+*>  X  [o,  »)  S- 

desired  properties.  | | | 


for  which  f  0  .  Then  t|>  :■ 
a  +  [-R,0]  with  a  M+R-i  . 
and  so  ip  ^  (  *— a )  has  the 


It  follows  from  Lemma  4  that  there  is  a  finite  collection  B  of  functions  supported 
on  [0,M]  whose  restrictions  to  [0,1]  form  a  basis  for  (S">|[0,  .  We  use  the 

abbreviation 

bj  ;■  b(»  -  J)  ,  b  €  B  ,  j  €  Z  . 

The  peeling  off  argument  of  Lemma  3  shows  that  the  functions 

bj(  b  e  B  ,  j-0,1,2,  ... 

form  a  basis  for  S"  .  But,  in  order  to  insure  local  linear  independence  of  this  basis,  we 
need  to  choose  B  with  some  care.  Precisely,  we  choose  B  so  that 

length  (B)  :■  length(b) 

is  minimal.  Here,  we  mean  by  length(b)  the  smallest  integer  r  for  which  supp  b  C 
[i,i+rj  for  some  integer  i  . 


Lena  5.  For  any  Interval  J  :»  [r,r+R]  ,  with  r  an  arbitrary  integer  and  R  as 
before,  the  collection 


*j  :=*  {bj  :  supp  bj  f]  J  f  I?  } 
is  linearly  Independent  over  J  . 

Proof.  Suppose  that  we  have  f  s“  c^ip  -  0  on  J  .  Then,  by  Lemma  3,  we  can 

J 

find  some  g  €  S”  so  that  f,  :*  f  -  g(*  -  (r+R>)  vanishes  on  [r,®)  .  If  not  all  the 
coefficients  c ^  were  zero,  then  we  could  choose  some  tp  :»  bj  with  c^  =*:  /  0  and 

so  that  j  is  as  small  as  possible.  Since  B  is  linearly  independent  over  [0,1]  , 

(f,)|  j+1j  f  0  ,  hence  we  must  have  J  <  r  .  Since  f^  vanishes  outside  [j,r]  , 

while  supp  bj  intersects  [j,j+1]  as  well  as  J  ,  this  would  imply  that 


lengthffy)  <  r-j  <  length(bj) 

while  f2  :«  fyl'+j)  €  S"  and,  on  [0,1],  is  in  the  span  of  B  with  the  coefficient  of 
b  not  zero.  But  then  B'  :=  (B\{b})U{f2)  would  also  supply  a  basis  and  length  (B ' )  < 
length(B)  ,  contradicting  our  choice  of  a  basis  of  minimal  length.  | | | 


Suppose  now  that  f  =  £  c^  ip  with  the  sum  over  all  ip  =  b^  with  j  €  Z  and  b  f  B  . 
Then  the  local  linear  independence  just  proved  allows  the  conclusion  that 
(2-2)  |cb.|  <  const  'f|[j(j+R)'„ 

for  some  constant  independent  of  b  and  j  . 

Since  S |  ^_R  qj  is  finite  dimensional,  we  can  choose  a  finite  dimensional  subspace 
T  of  S  which  is  carried  faithfully  onto  S| j_R  qj  by  the  restriction  map 
f  | — >  f  | [_R  o]  •  For  s  €  S  ,  let  sT  be  the  unique  element  of  T  which  agrees  with 
s  on  [-R,0]  .  Then 

(s-aT^  X[o,  €  S 

by  Lemma  3,  hence  on  [0 ,<*>) , 


(2.3) 


*  “  ®T  +  £j>0  ^bSB  cbj  bj 


for  suitable  coefficients  cfc  .  Since  s |  [ o , <*»)  is  inf*nite  dimensional  by  (2.1.v),  this 
show3  that  B  ?  0  and  that  S~  is  also  infinite  dimensional.  Further,  since  the  linear 

2 

space  of  all  linear  maps  from  T  to  S|  j_R  gj  has  dimension  (dim  T)  ,  the  maps 

T  _>  S|[-R,0):  f  l_>  f<*  ~  a)|(-R,0]  '  a  f  A 

must  be  linearly  dependent  whenever  A  C  Z  and  ] A |  >  (dim  T)2  .  This  together  with  Lemma 

3  implies  that,  for  any  A  C  Z  with  more  than  (dim  T)^  elements,  there  exist 


(2.4) 


1  ,  and 

X[0,«)  £a£A  caf(’  -  a)  €  s’  for  a11  f  €  T  • 


Proof  of  Theorem  1.  Let  g  be  the  continuous  piecewise  linear  function  with  integer 
breakpoints  which  has  the  value  1  at  x  =  1,  2,  3  and  vanishes  at  all  other  integers. 
Since  g  C  Cq(R)  ,  there  is  sb  5  Sh  such  that  Sy,  converges  to  g  uniformly  on  R  . 
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From  now  on,  w«  consider  only  h  €  1/M  .  We  write  sb  t®  o^*h  ®n<*  9h  °h?  •  ^*n 

«gh  -  »h«  — >  0  as  h  — >  0  .  Since  sh  €  S  ,  we  have  sh  *  th  +  sb  on  [0,®)  ,  with  th 

€  T  and  sJJ  6  s"  •  Now  choose  A  C  8H/h  as  above.  Then  the  functions  gh(  •  -  a)  ,  a€A  , 

are  supported  in  [0,®)  and  have  disjoint  supports.  We  define  Gh  :•  IagA  ca  g^(*-a)  , 

with  the  ca  as  given  by  (2.4).  Then  fh  t®  ca  sb(*-a)  is  in  S'  ,  and 

#f.  -  G.  I  — >  0  as  h  — >  0  . 
h  h 

Now,  because  of  (2.4),  Gj,  is  identically  one  on  an  interval  of  the  form  [mh-1/h, 
mh+1/hj  .  Hence,  with  Fh  i«  fj,(*+mh)  ,  we  have 


,Fh  •  1,M/h,l/h]  0  “8  h“>  ° 


Since  fh  €  S',  we  can  write 

In  view  of  (2.2), 

(2.5) 


ph  Ej€s  ^b^B  cbj(h)  bj  * 


| cb ^  < h ) |  <  const  for  all  b,  j,  h  . 

We  can  take  therefore  a  subsequence  h’  of  h  so  that,  for  all  b  and  j  ,  Cj^(h') 

converges  as  h*  — >  0  to  some  cb^  .  It  follows  that 

Ej€*  ^  cbj  bj 

converges  uniformly  on  compact  seta  to  some  function  F  €  S  .  Since  we  also  have  that 
Fh  converges  to  F  uniformly  on  compact  sets,  we  must  have  F  -  1  .  Thus 

(Cj^bj )bgB( is  a  partition  of  unity  in  S  .  It  is  good  because  of  (2.5).  ||| 
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We  show  that  for  certain  translation  invariant  spaces  S,  a  necessary 
and  sufficient  condition  for  the  eventual  denseness  of  the  corresponding 
scaled  spaces  is  that  S  contain  a  stable  and  locally  supported 

partition  of  unity. 

These  results  have  been  motivated  by  recent  work  on  approximation  oy 
multivariate  piecewise  polynomials  on  regular  meshes. 


DO  /.ST,*  1473  EDITION  OF  1  NOV  «»  1$  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  f1Wi»n  Data  Entarad) 


I 


